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Entanglement from the Dynamics of an Ideal Bose Gas in a Lattice 
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We show how the remotest sites of a finite lattice can be entangled, with the amount of entan- 
glement exceeding that of a singlet, solely through the dynamics of an ideal Bose gas in a special 
initial state in the lattice. When additional occupation number measurements are made on the 
intermediate lattice sites, then the amount of entanglement and the length of the lattice separating 
the entangled sites can be significantly enhanced. The entanglement generated by this dynamical 
procedure is found to be higher than that for the ground state of an ideal Bose gas in the same 
lattice. A second dynamical evolution is shown to verify the existence of these entangled states, as 
well entangle qubits belonging to well separated quantum registers. 
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One of the aims in the field of quantum information is 
to set up entanglement between locations separated by 
some distance, and in general, greater the separation and 
more the amount of this entanglement, the better. While 
photons are best for long distance entanglement distribu- 
tion, for short distances (such as for linking quantum reg- 
isters) other alternatives are important ^ n this con- 
text, the dynamics of spin chains have been proposed for 
the distribution of entanglement over a distance of sev- 
eral lattice sites (For example, Refs. 0,130] to mention 
a very few). However, as the number of possible states 
of a spin in a spin chain is low, the amount of entangle- 
ment that can be dynamically generated and distributed 
through a single spin chain channel in a limited time is 
restricted. In this paper, as an alternative to spin chains, 
we suggest the use of the dynamics of an ideal gas of M 
bosons in a lattice to generate and distribute entangle- 
ment between its remotest sites. Note that our dynamics, 
when followed by certain occupation number measure- 
ments, will create high but "finite" dimensional entangled 
states (with entanglement ~ log 2 V^+Const.) which 
are qualitatively very different from the infinite dimen- 
sional Gaussian entangled states which can be generated 
dynamically through harmonic oscillator chains [f| . 

Another motivation for the current work originates 
from the literature on entangling Bose-Einstein conden- 
sates (BECs) of gaseous atoms/molecules in distinct 
traps ( Refs. 0, 0, IE 13 113 to mention a very few) where 
usually small lattices or continuous variable entangle- 
ment are considered. Can we use lattice dynamics to 
create entanglement between traps separated by several 
intervening lattice sites'? Here we show how to accom- 
plish this without either the physical movement of traps 
or any local modulation of the lattice parameters. 

We consider a one dimensional lattice of N sites, where 
the aim is to establish a significant amount of entangle- 
ment between sites 1 and N. We choose N to be odd and 
initially place M bosons in the ^jpth site of the lattice 
and keep all the other sites empty. Physically this cor- 
responds to a Fock state \M) on the ^j^th site of the 
lattice and a vacuum state in all the other sites. In terms 
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FIG. 1: Our setup of creating entanglement between sites 1 
and N of a ID lattice. One simply starts with M bosons in 
the (N+ 1) /2th lattice sites and allows dynamical evolution of 
the system to create entanglement between the sites 1 and N. 
Alice and Bob who have access to the sites 1 and N can use 
the entanglement for quantum communications of for linking 
distinct quantum registers. 



of boson creation operators which create a boson in 
the jth lattice site the initial state is thus 
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This special initial state greatly simplifies the calcula- 
tions and can be regarded as a generalization of a single 
spin flip at the midpoint of a spin chain, which has been 
studied for entanglement generation We assume the 
bosons to be essentially non-interacting during the time- 
scale of our scheme i.e., their collection is an ideal Bose 
gas. Then the Hamiltonian of the system in the lattice is 
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H = J ^(atoj- + i + Oja) +1 ). 
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The state of each boson then evolves independently in the 
lattice {i.e., each boson evolves as if it was hopping alone 
in an otherwise empty lattice) as a «+i — > fj(t) a \ 
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where fj(t) is the amplitude of the transfer of a single 
boson from the ^^th site to the jth site in time t. Thus 
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the state of the M boson system at time t is 



(3) 



The evolution amplitudes /j (t) are identical to that of a 
XY spin chain in the single excitation sector, and is given 
HQ by 



— > J sin — sin— le "+ 1 . (4) 
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Eqs.© and J3J give the complete time evolution of the 
M boson state analytically. 

Our task is now to calculate how much entanglement 
exists between sites 1 and N in the state of Eq.Q as 
a function of time and find a time at which this entan- 
glement is large. We will thus have to calculate the re- 
duced density matrix of the states of sites 1 and N by 
tracing out the state of all the other sites. To accom- 
plish that we adopt a strategy from Ref . 01 > which eval- 
uates the entanglement between two regions of a Bose- 
Einstein condensate in its ground state in a single trap. 
The strategy is to define new creation/annihilation oper- 
ators by combining the operators as = fi(t)(a\ + 

ajv)/V2FiWF and L+ = fj (t)o}/ y/l - 2\f } (i)| 2 , 

which are valid bosonic creation operators satisfying 
[E, E^} = 1, [L, L^] = 1. Noting that the symmetry of 
our problem implies fi(t) = /at(£), we expand the expres- 
sion of \^(t)) in Eq.© in terms of the above operators 
to get 
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|tt(t)) - £ A /^(v/2|/ 1 (i)|n A /l-2|/ 1 (i)| 2 )' 

r=0 

X |^r)liv|0r)2...JV, (5) 



where we have substituted \ip. 
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|^r/2...jv = y) ®'' Noting that the sc ^ °^ states 

{\4>r)2...N} represents an orthonormal set in the space of 
states of the sites 2 to N — 1, we have the the reduced 
density matrix of sites 1 and N to be 
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where P r (t) = M C r (2\f 1 (t)\ 2 r(l-2\f 1 (t)\ 2 ) M - r . We 
can write \i/) r )iN i n terms of the occupation numbers as 
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Note that the only time dependence of the state /o(t)uv 
stems from fi(t), maximizing which over a long period of 
time is a pretty good strategy for obtaining a time th such 



that p(th)iN is highly entangled. The maximization en- 
sures that the proportion of the state |^jw)iiVj which has 
the most entanglement among the set of states {|Vv)iat}, 
is the highest possible in p(t h )i N . Ideally, the lattice 
dynamics should be frozen at t^, say by globally rais- 
ing the barriers between all wells of the lattice, so that 
Alice and Bob can utilize p{th)\N for quantum commu- 
nications or linking quantum registers. For the smallest 
non-trivial lattice [N — 3) we know that fx(t) = 1/V2 
at t = 7r/2J\/2 from the XY model [llj. For this case, 
the state \iPm) 13 is generated between sites 1 and 3 at 
t = 7r/2JV2 whose entanglement can be made to grow 
without limits by increasing M. For this special case, the 
advantage over spin-1/2 chains is most evident, where 
only the case of M = 1 can be realized (with a single flip 
in the middle of a chain of three spin-1/2 systems). 

In general (for N > 3), the state p(t)\N is a mixed state 
of a (M + 1) x (M + 1) dimensional system, and the only 
readily computable measure of its entanglement is the 
logarithmic negativity E n [l2T | , which bounds the amount 
of pure state entanglement extractable by local actions 
from the state p{t)i^. It is the standard measure used 
when high dimensional mixed entangled states arise 
The E n of p(t)w for an appropriately chosen time t = th 
(see above) are plotted in FigEl as a function of N for 
different values of M. This figure clearly shows that for TV 
as high as 21, one can generate more entanglement than 
that of a singlet (E n = 1 for a singlet). For such modest 
lengths, thus one generates more entanglement between 
the ends of a lattice by using a M > 1 boson gas than 
ever possible with a spin-1/2 chain, which is the M=l 
case, also plotted in FigEl The advantages of increasing 
M diminish, though, as one increases N. For N = 51, we 
see that though there is some advantage of high boson 
number (M = 50) over the spin chain (M — 1) case, 
this advantage does not increase by increasing M (for 
example, the M = 5 and 50 plots are nearly coincident). 

Next, we slightly modify our scheme and after the dy- 
namical evolution till th, we measure the total number of 
bosons in the intermediate sites (sites 2 to AT — 1). With 
probability P r , we will find this number to be M — r, and 
when we do so, we will generate the state \ipr}iN between 
sites 1 and N. Note that \tp r )iN is created between sites 
1 and N whenever a "total" of M — r bosons is found 
in the remaining sites irrespective of the distribution of 
these M — r bosons among the sites. The amount of 
entanglement in the pure state \ip r )iN is given by its 
von Neumann entropy of entanglement E v (\ip r )iN) = 
-Tr(pilog 2 pi), where pi = Tr N (\ip r )(tp r \ 1N ), which 
equals the quantity of entanglement that can be obtained 
as singlets (the most useful form, say for their use as a 
resource for perfect teleportation of qubits from Alice to 
Bob) from the state by local actions and classical commu- 
nications alone. Thus the average von Neumann entropy 
of entanglement (E v ) = ^2 P r E v (\tjj r )iN) over all pos- 
sible measurement outcomes is operationally the most 
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useful measure of the entanglement between sites 1 and 
N in our modified scheme. This quantity has been plot- 
ted in Fig|3|for various M for large lengths of lattice up 
to N = 1001. We find that for M = 1000, entanglement 
nearly equal to that of 4 singlets is generated across a dis- 
tance of 1001 lattice sites, and this is more than 70 times 
the amount possible with a spin-1/2 chain (M = 1) of 
same length. For high N and M, we can represent P r by 
a Poisson distribution and know that /i(t) ~ 1.7/N 1 ' 3 
at t ~ (N + 0.817V 1/3 )/4J HEl Using t hese, w e have 
the analytic expression (E v ) ~ \og 2 {1.7y / MTre/N 1 / 3 }, 
whose fit with data gets better as we proceed from the 
M = 10 to the M = 1000 plot at high N. For every or- 
der of magnitude increase of M we thus expect to gain an 
entanglement equalling that of log 2 \/l0 = 1.66 singlets, 
which is confirmed by the differences in (E v ) between the 
M = 100 and 1000 plots at high N (FigE|). So for ob- 
taining entanglement exceeding that of 10 singlets across 
a distance of ~ 1000 lattice sites we require M ~ 10 7 . 

Now we proceed to compare the amount of entangle- 
ment between sites 1 and N generated by the above 
schemes with that obtainable from the ground state of an 
ideal Bose gas in the same lattice (to check whether we 
have gained from dynamics). Ground state entanglement 
between distinct regions of a Bose-Einstein condensate 
in a single trap has already been investigated [Hil Il4| . 
In our lattice setting, the ground state of H is simply 
(7^ J2f=i a]) M \Q}/VM. It thus suffices to replace /i(t) 

in the expression of P r (t) with l/y/N to change from the 
dynamical state to the ground state. E n between sites 
1 and N for the ground state is plotted using asterisks 
in FigH as a function of N for M = 50. We find that 
the entanglement of sites 1 and N is nearly vanishing for 
N > 25, and thus there is a marked advantage of using 
dynamics as opposed to the ground state. One may also 
however, measure the occupation numbers of sites 2 to 
N of the ground state to obtain an average entanglement 
(E v ) ~ log 2 \J Mne/N between sites 1 and N (for high 
M and N). Even this is lower than that of our second 
scheme by log 2 {1.7N 1 / 6 }. 

We now proceed to discuss a method to verify the 
pure entangled states \ip r )iN created through our sec- 
ond scheme. One can easily verify the number corre- 
lations (the fact that if site 1 has k bosons, then the 
site N has r — k bosons) in \ip r )iN by occupation num- 
ber measurements. Then the only task remaining is to 
verify the coherence between the terms |fc)i|r — k) ^ in 
|Vv)iiv- If we a Pply a Hamiltonian H v = na\ai for a 
fixed time t v to site 1 then the above coherence leads to 
|Vv)ijv becoming \ip' r )iN = {e l8 a\ +aj v ) r |0)/V2 r r! where 
= Kt v . The ideal Bose gas is then allowed to evolve 
again in the lattice (assuming the lattice sites 2 to TV — 1 
have been emptied due to or after the earlier occupation 
number measurements), which results in its state becom- 
ing |*' (<)>!..* = (Ef=i{e J V-W+.9^«)}«])''|0)/VrH 
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FIG. 2: Entanglement (log negativity) of sites 1 and JV from 
dynamics at an optimal time for various iV and M, and for 
the ground state. 



where gij are amplitudes for a boson going from the Zth 
site to the jth site. From symmetry and the formulae for 
evolution in XY chain, we have q x n+i (t) = q n n+i (t) = 

fi (t) . Thus the probability of finding the site oc- 
cupied can be varied from 1 — (1 — 2\fi(t)\ 2 ) r for 6 = 0, 
to for 8 = 7r (as the term a/ N+1 in | \T/ (t))i..jv vanishes 
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for = 7r). This variation with 9, whose range increases 
with increasing r and which can be further increased by 
maximizing enables us to verify the coherence be- 

tween the terms |fc)i|r — fe)iv in \ip r )iN- One way to verify 
the mixed state p(t)iN produced by our first scheme will 
be to go through our second scheme and check that pure 
states \4> r )iN are produced with probabilities P r . 

Next, we provide an example of linking qubits A and 
B of distant quantum registers by making them in- 
teract with sites 1 and N respectively for a time t q , 
when these sites are in the state \ip r )iN- The initial 
state of the qubits is {\0) A + |1)a) ® (|0) B + |l)fl)/2 
and the interaction Hamiltonian is H q = g(afa\a\ + 
afa^a^) in which and af are Pauli operators of 
the qubits. For gt q = n, the state of the system is 

{(\0)A\0)B + (~iy\l)A\l)B)®(a\+a i l ) r \0) + (\0)A\l)B + 
(-l) r |l).4|0) B ) ® (a\ - aj v ) r |0)}/\/2 T +2H. As in the 
previous paragraph, now the state of the ideal Bose 
gas is again allowed to evolve in the lattice and after 
some time the presence of any boson the (N + l)/2th 
site is measured. If this site is found occupied, then, 
remembering the logic of the last paragraph, the gas 
must have been in the state (a\ + aj v ) r |0)/V2 r r!, which 
projects the qubits to the maximally entangled state 
|0) J 4|0) B + (-l) r |l) J 4|l)i3- The probability for this is 1/2. 
Otherwise (also with probability 1/2), a state whose fi- 
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FIG. 3: Average entanglement (von Neumann entropy) of 
sites 1 and N from dynamics followed by measurements on the 
intermediate lattice sites for various N and M. The plot with 
asterisks denotes the same for the ground state for M = 1000. 



delity with |0)A|l)s + (-l) r |l).4|0)B is 1 - (1 - 2 1 A (t) | 2 
is obtained (this fidelity can be maximized by maximizing 
|/i(t)|) and gets better with increasing r). This scheme 
for entangling qubits of distinct registers is just an ex- 
ample. The values of (E v ) found earlier imply that in 
principle one should be able to extract many more sin- 
glets from \ip r )iN by local actions alone. 

BECs of dilute atomic gases in optical lattices form a 
test ground for our protocols. Fock states can be pre- 
pared by Mott transitions |l5j, the interactions needed 
for that can be switched off by a Feshbach resonance 
[lif to obtain an ideal Bose gas (otherwise, one sim- 
ply has to go between regimes were the on-site repulsion 
U » J/M 2 toll « J/M 2 N by the global modulation 
of lattice potentials). Accurate measurement of the total 
number of bosons in sites 2 to N ~ 1 is possible either by 
using metastable atoms |l7j or potentially through a sec- 
ond Mott transition involving these sites only (note that 
individual site occupation numbers are not required) . For 
the verification part, only whether the (N + l)/2th site 
is occupied or not need to be ascertained, and atomic 
fluorescence in external fields can potentially be used. 
By placing all atoms in a known magnetic state when 
required (say by a laser), each site can be imparted a 
magnetic moment proportional to its occupation number. 
This, and a local magnetic field at site 1 can be used to 
realize H v . This magnetic moment also enables magnetic 
moment based register qubits (atomic or solid state) to 
interact with sites 1 and N through H q . If imparting 
the atoms with a magnetic moment automatically make 
U >> J/M 2 then we have to ensure that Ut v and U t q are 
integral multiples of 2tt. An alternative implementation 



is nano-oscillators arrays |5| when resonantly coupled to 
each other. They can be coupled to Cooper-pair box 
qubits to both create and measure Fock states [lSj . and 
to Josephson qubits through a Jaynes-Cummings model 
whose off resonant limit can implement H q . Coupled 
cavities in photonic crystals, where Fock state prepara- 
tion and measurements could be performed with dopant 
atoms is another potential implementation |2(lj . 
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